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Low energy excitations in carbon nanotubes can be de- 
scribed by an effective field theory of two components spinor. 
It is pointed out that the chiral anomaly in f+1 dimensions 
should be observed in a metallic toroidal carbon nanotube on 
a planar geometry with varying magnetic field. We propose 
an experimental setup for studying this quantum effect. We 
also analyze the vacuum structure of the metallic toroidal car- 
bon nanotube including the Coulomb interactions and discuss 
some effects of external charges on the vacuum. 



I. INTRODUCTION 

In recent years carbon nanotubes(CNTs) |Q] have attracted 
much attention from various points of view. Especially their 
unique mechanical and electrical properties have stimulated 
many people's interest in the analysis of CNTs |^-^. They 
have exceptional strength and stability, and they can exhibit 
either metallic or semiconducting depending on the diameter 
and helicity |5|Q. Because of their small size, properties of 
CNTs should be governed by the law of quantum mechanics. 
Therefore it is quite important to understand the quantum be- 
havior of the electrons on CNTs. The bulk electric properties 
of (single- walled) CNTs are relatively simple, but the behavior 
of electrons at the end of a tube(cap) or metal-CNT junction 
is complicated and its understanding is necessary for build- 
ing actual electrical devices. On the other hands, toroidal 
CNTs (FuUerence 'Crop Circles' Q) are clearly simple be- 
cause of their no-boundary shape and they can also be either 
metal or semiconducting properties( hereafter we use 'torus' 
or 'nanotorus' instead of 'toroidal carbon nanotube' for sim- 
plicity). Even in the torus case, quite important effect: 'chiral 
anomaly' which is of essentially quantum nature, might 
occur. 

Low energy excitations on CNTs at half filling move along 
the tubule axis because the circumference degree of freedom 
(an excitation in the compactified direction) is frozen by a 
wide energy gap. Hence this system can be described as a 
1-1-1 dimensional system. Furthermore in the case of metallic 
CNTs, the system describing small fluctuations around the 
Fermi point is equivalent to the massless fermion in 1+1 di- 
mensions. If we include gauge field, this situation can be mod- 
eled by the quantum field theory of massless fermion which 
couples to the gauge field through minimal coupling. This 
model realizes the chiral anomaly phenomenon |9|-[ll[. 

The chiral anomaly is one of the most interesting phenom- 
ena in quantum field theory and has had an appreciable influ- 
ence on the modern development of high energy physics |l^] 
and of condensed matter physics . The effect of the chiral 
anomaly on the electrons in a nanotorus appears directly as 
a current flow. On the other hand, it is known in solid state 
physics that a one-dimensional metallic ring shows the 'persis- 
tent current [Q' in an appropriate experimental setting. The 



current originating from the chiral anomaly shows the same 
magnetic field dependence to the persistent current. There- 
fore, the persistent current is closely connected with the chiral 
anomaly in 1-1-1 dimensions. The chiral anomaly provides a 
deeper understanding for the persistent current as shown in 
the present paper. 

In this paper, we examine the anomaly effect in a metallic 
nanotorus and discuss how such an effect can be observed 
experimentally. We also clarify the vacuum structure of the 
model regarding the gauge field as a classical field and discuss 
some effects of external charges situated on the metallic torus. 

The organization of this paper is as follows. After reviewing 
quantum physics of CNTs 1 15 , we study the case of nanotorus 
in section O. In section [II we point out that low energy 
excitations on a metallic torus at half filling can be modeled 
by a quantum field theory of massless fermions with gauge 
field and construct the Hamiltonian of this system. In section 
[[V| we discuss the chiral anomaly and show how such an effect 
can be observed experimentally. We examine the Hamiltonian 
including the Coulomb interactions and analyze an effect of 
the Coulomb interactions on the chiral anomaly in section 
We discuss some effects of external charges in s ection VI 
Conclusion and discussion are given in section VH. 



II. CARBON NANOTORUS 

A carbon nanotube can be thought of as a layer of graphite 
sheet folded-up into a cylinder. A Graphite sheet consists 
of many hexagons whose vertices are occupied by the car- 
bon atoms and each carbon supplies one conducting electron 
which determines the electric properties of the graphite sheet. 
The lattice structure of a two-dimensional graphite sheet is 
shown in Fig.|l| 
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Lattice structure of a two-dimensional graphite 
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It is obvious, however, that this picture of a CNT as a 
graphite sheet rolled up to form a compact cylinder is some- 
what oversimplified. We need to be careful of the follow- 
ing facts. First, a conducting electron makes the tt orbitals 
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whose wave function extends into the Z-direction: perpendic- 
ular to the graphite sheet. Hence, in the case of muhi- walled 
CNTs(MWCNTs), the wave functions which belong to differ- 
ent layers may interfere so that there is a chance that some 
electrical properties will alter as compared with SWC- 
NTs. Second, we should care for the curvature of a cylinder. 
This cause the mixing of a and tt orbital so that the band 
structure might change. In the present paper, we concen- 
trate single-walled CNTs with a large diameter. Therefore 
the effects of interlayer interactions and curvature can safely 
be neglected [ p^ . 

There are two symmetry translation vectors on the planar 
honeycomb lattice, 



%/3 , 3 
— ae^, + -^aey 



(2.1) 



Here a denotes the length of the nearest carbon vertex and 
is given by 0.142 nm. and ey are unit vectors which are 
orthogonal to each other [e^ ■ ey = 0). If we neglect the spin 
degrees of freedom, because of these symmetry translations, 
the Hilbert space is spanned by the following two Bloch basis 
vectors, 



lvl>^) = ^e"='--al|0), 
l*^) = ^e""--al|0). 



(2.2) 
(2.3) 



where the black(») and blank(o) indices are indicated in Fig.^. 
Ti labels the vector pointing each site i, and ai,aj are canon- 
ically annihilation-creation operators of the electrons of site i 
and j that satisfy 



(2.4) 



a cylinder by imposing a periodic boundary condition to the 
state vector. In general, we may consider the following bound- 
ary condition. 



G(iVri + MTa 



G denotes a symmetry translation operator. From which we 
obtain 



^ {2N + M) aK + '^Maky 



■ 27rn, 



(2.9) 



where n is an integer. Next, in oder to make a torus, we 
compactify the tube into a torus by imposing a boundary 
condition to the tubule axis direction. For example, we may 
consider a 'zigzag torus' which has the following boundary 
conditions. 



G{M{2T2 -Ti))!*'') 



(2.10) 



The former condition makes a sheet into a zigzag tube, and 
the latter forces the tube into a zigzag torus. 

It is clear that there are many possibilities for the shape 
of torus and each shape has its own boundary condition. So, 



some of them have different properties from the one (2.10). 
Especially we can image a torus in which some twist exists 
along the tubule axis direction |l^. This system have the 
following boundary conditions. 



G(Afri)i*''') = I**), 

G{M{2T2 - Ti))!*'-') = G(iVri)|*'''), 



(2.11) 



where A'' is decided by the twist at the junction of tube end, 
see Figj^. These boundary conditions yield the discrete wave 
vectors 



We construct a state vector which is an eigenvector of these 
symmetry translations as follows: 



j^fc^ = c.*i*i) + 6-01*0). 



(2.5) 



In order to define the unit cell of wave vector fc, we act the 
symmetry translation operators on the state vector and obtain 
the Brillouin zone 

TT TT \/3 ,3 /„ „\ 

1= < akx < —TV < ——akx + -aky < tt, (2.6) 

where k^ — k ■ e^. 

The mapping of the graphite sheet onto a cylindrical sur- 
face is specified by a wrapping vector: 



C = NT1+MT2, 



(2.7) 



which defines the relative location of the two sites. The pair 
of indices (A*', Af) describes how the sheet is wrapped to form 
the cylinder, and determines its electrical properties [^,^. It 
can be divided CNTs into three categories depending on the 
pair of integers (A'^, M). A tube is called 'zigzag' type if M = 
and 'armchair' in the case M — —2N. All other tubes are of 
the 'chiral' type. 

To study the electronic properties of a nanotorus quan- 
tum mechanically, first we compactify the graphite sheet into 



akx 



2-K n 



""^^ 3M 



m + 



(2.12) 
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FIG. 2. Twisted zigzag torus: We impose the periodic 
boundary condition in the X-direction and attack same char- 
acter('B') in the Y-direction to obtain a twisted nanotorus. 
There are A'' hexagons between the point A and B. 
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So far we have constructed the Hilbert space of the con- 
ducting electrons. The Hilbert space is spann ed by the Bloch 
basis vectors with the discrete wave vectors (2.12). Now, we 
consider the Hamiltonian which governs the time evolution 
of a state vector. Each carbon atom has an electron which 
makes vr-orbital. The electron can transfer from any site to 
the nearest three sites through the quantum mechanical tun- 
neling or thermal hopping in finite temperature. Therefore 
there is some probability amplitude for this process. In this 
case, the tight-binding Hamiltonian is most suitable: 



(2.13) 



where the sum is over pairs of nearest-neighbors carbon 
atoms i,j on the lattice. 7 is the transition amplitude from 
one site to the nearest sites and -Eo is the one from a site 
to the same site. This parameter Eo only fixes the origin of 
the energy and therefore irrelevant. It is an easy task to find 
the energy eigenstates and eigenvalues of this Hamiltonian. 
In the matrix representation, the energy eigenvalue equation 
reads 



Eo 



Eo 



where 



1*0) = 



Ek 



c: 



(2.14) 



(2.15) 



and the vector Ui is a triad of vectors pointing respectively 
in the direction of the nearest neighbors of a black point, see 
Fig.0. The energy eigenvalues and eigenvectors are as follows, 



Ek = Eo± A(fc) 

c':\ 1 



^A(fc) 



Eiku 

±A(fc) 



(2.16) 
(2.17) 



where 



3 \/3 

A(fc) = 7\/ 1 -I- 4 cos akx cos —aky + 4 cos^ ak^. 

V 2 2 2 

The structure of this energy band has striking properties when 
considered at half filling. This is the situation which is physi- 
cally interesting. Since each level of the band may accommo- 
date two states due to the spin degeneracy, the Fermi level 
turns out to be at midpoint of the band(ii'fc = -Eo). Fermi 
points in the Brillouin zone are located at 



cifci,2 = [akx, aky) 



,^ 2-K 2tt. 



(2.18) 



Hence, if A'^ in eq.(2.12) is a multiple of 3, then the zigzag 
torus shows metallic properties. 

In order to understand the electric properties we should 
take into account a small perturbation around the Fermi 
point. So we take k = ki + Sk as a. small fiuctuation | ^ . 
Perturbation around the point k2 is same as around the point 
k = k\. So we may only consider one of the pairs. In this 
case the Hamiltonian which describes the system is given by 



-70- • 5k ■ 



-7a 



Sk, 





+ iSky 



5kx 



iSky 







(2.19) 



This implies that the low energy excitations of metallic CNTs 
at half filling are described by an effective theory of two di- 
mensional spinor obeying the Weyl equation. In coordinate 
space, we find 



U 



pert 



37a 



■a ■ p, 



(2.20) 



where p is the momentum operator, p — —ihV and ai are 
the Pauli matrices. It is also convenient to use a parameter 
f3 = — and T = f3t. In this case the Schrodinger equation 
becomes l20|l 



d 



(2.21) 



In the following section, we consider metallic and semi- 
conducting zigzag torus that have small N and large M 
values(^ > 10^). In this case, transitions between difi'er- 
ent kx are very small because of their costed energy (~ ■'^) as 
compared to that of ky:{^ J^). Therefore, the only surviv- 
ing degree is a motion in the tubule axis direction, i.e. this 
system is a 1-1-1 dimensional system effectively. 



III. EFFECTIVE FIELD THEORY OF CARBON 
NANOTUBE 

In this section we would like to focus on the zigzag torus 
which has the boundary conditions ( ^.ICj ) and construct an ef- 
fective field theory describing the low ener gy exc itations in the 
torus. More general boundary conditions(2.11 ) are discussed 
in section IV , The zigzag torus can exhibit either metallic or 
semi-conducting depending on the value of A''. If A*' is a mul- 
tiple of 3 then the torus shows metallic properties. In order to 
analyze the semi-conducting case equally, we set N = ±3n-|-6, 
where n is a positive integer and b G {0, ±1, ±2}. To examine 
the low energy excitations, we should consider the following 
wave vectors and energy: 



akx = ± 
Ek = Ek 



2-K 



1 



27r 

, aky — + adky, 

3\/3 1±^ 3 



3ti 

2-K 1 



2tt 



— ,T— +aSky 



(3.1) 
(3.2) 



Considering the small perturbation {aSky < ^) and a large 
diameter (n > 3), the excitation energy can be approximated 
by 



Ekr^±VFy^AP+P^, (3.3) 
where «f(= |/3|) is the Fermi velocity, p = hSky and 
2nh \b\ 



M ■ 



(3.4) 



Consequently, we have obtained a linear dispersion relation 
for metallic case(6 = 0). This is to be contrasted with the 
dispersion relation for semi-conducting case(6 7^ 0). 
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Here we comment on the excitations between different k^. 
In order to neglect these excitations, we have to restrict the 
energy to the following region: 



Ek < Esh^ = 



7- 



(3.5) 



At room temperature, thermal excitation between different 
\5kx\ cannot occur by the Boltzmann suppression factor, 

^Skj^ 200 

e ET-(r.^ e~~f"), so that the effects of different \5kx\ can 
be ignored. Accordingly, the Hamiltonian which describes 
the low energy excitations near the Fermi point is given by 



37a / ±M- 
2^1 ±A/ + ip 



ip 



(3.6) 



The Coulomb interaction term is 

^2 >■ J\x)j\y) 



-He = 



Stt 



\x ■ 



■y\ 



dy. 



(3.13) 



where ej''{x) stands for the charge density. The electric cur- 
rent is given by ej^ = e#7^^' where e is the electron charge. 

It should be note that besides the Coulomb interaction, 
backscattering and umklapp process may come in the dynam- 
ics of electrons as is shown in reference p3| . We neglect these 
interactions providing that their coupling are weak enough. 



IV. CHIRAL ANOMALY IN A METALLIC 
NANOTORUS 



The sign (±) in front of M comes from the ki and k2 points. 
The sign can be removed by an appropriate unitary trans- 
formation of the state vector. Therefore we may choose the 
minus sign without loss of generality. Using the unitary trans- 
formation, U = e"'^"^"", the Schrodinger equation becomes 



P 

-M 



-M 
-P 



(3.7) 



One can obtain the quantum field theory by promoting 
the wave function(7/)) to the field operator(\[') obeying the 
anti-commutation relation. Because the Schrodinger equa- 



tion (3/7) is the Dirac equation in two dimensions, it is ap- 
propriate to adopt the following Lagrangian density. Addition 
of the electro-magnetic interaction according to the minimal 
coupling gives 



C = - -F^'F,,, - t (^ -h M) 



(3.8) 



where D is the covariant derivative and is the field 

strength. 



fi=0,l 



(3.9) 



The gauge fields(A^) live in four dimensional space-time so 
that summention indices run from to 3 in the gauge kinetic 
term. We adopt the following notation; 

g'"' = diag{l, -1,-1,-1}, 7"=^'j' > 



1 



.5 .01 



-10/'^ ^ \ -1 



1 



(3.10) 



We quantize the fermion field in each configuration of 
the gauge field choosing Weyl gauge condition {Ao = 0). 
Fermionic part of the Hamiltonian density is given by 



Ti, = Hf + TLc- 



(3.11) 



The total Hamiltonian density consists of the kinetic term 
and the Coulomb term. The kinetic term is given by 



nF= -a^^hp-^ 



ihdi - ^Ai 

^ c ^ 

-M 



-M 

-{ihdi - f Ai) 



(3.12) 



We have obtained the Hamiltonian density which describes 
the low energy excitations in the zigzag torus. The Hamilto- 
nian consists of two parts, one is the kinetic term and the 
other is the Coulomb interaction. In this section we dis- 
cuss the quantum mechanical vacuum structure of the kinetic 
Hamiltonian: Hf{= ^Hf)- Effects of the Coulomb interac- 
tion will be considered in latter sections. Hereafter we focus 
on the metallic case and set = 1 in this section for sim- 
plicity. In the metallic case: 6 = 0, the energy eigenvectors 
are given by 



tlFlpn 



hF^^n, 



— e Jo 



(4.1) 



where e„ is the energy eigenvalues and L is the circumferen- 
tial length of the zigzag torus: L — 3a\M\. We expand the 
fermion field using the left and right moving waves as 



■f{x,T) = 'fL{x,T) + -il'R(x,T) 



E 



a„il}„{x)e ' « 



-I- h„ip„{x)e 



(4.2) 



where a„, fe„ are independent fermionic annihilation operators 
satisfying the anti-commutators: 



(4.3) 



All the other anticommutator vanish. 

In order to get the energy spectrum e„, we have to i mpose 
a boundary condition to the eigenfunctions. We take (2.11) 
as a general boundary condition. The boundary condition of 
the zigzag torus in the Y-direction becomes the following: 



G{M{2T2 - Ti))!*'') = 6=^"^^!*''), 



(4.4) 



where plus(minus) in exponent has its origin in the Fermi 
point, ^1(^2). Hence we should impose the following bound- 
ary conditions on the fermion energy eigenfunctions. 



i,4x + L)^e^'-''i,n{x), 
then the energy eigenvalues are given by 
TV 



2Tv{n ± 



e 
he 



Aidy 



(4.5) 



(4.6) 
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Here, the gauge field Ai is experimentally controllable by the 
following experimental setup, see Fig.|[ On a planar geometry 
we set a nanotorus and put some magnetic field inside the 
torus perpendicular to the plane. In this case the gauge field 
that expresses this magnetic field is given by, in the vector 



notation, A 



N<t> 



p-\/9 where 6 is an angle of two points on 



the torus. Therefore we get a component, 
, N^D 

Ai = — z — , 



(4.7) 



where 4>d = 22^, This vector potential expresses the N flux 
inside the torus and by tuning the magnetic field, N can be 
taken as a real number. 

, Magnetic Flux 




■ Toroidal CNT 

FIG. 3. A toroidal carbon nanotube on a planar geometry 
with magnetic field 

The Lagrangian density ( ^.^ ) with massless fermion has 
two classical conserved currents 



(4.8) 
(4.9) 



where e''" is an antisymmetric tensor: e"^ = 1. Therefore, 
the following two charges conserve in the time evolution of 
the system at the classical level. 



j'^{x)dx, 



= <h Jl{x)dx. 



(4.10) 



(4.11) 



J'^ (electric current) conservation((?^ J'' = 0) is due to the 
gauge symmetry and (chiral current) conservation(9^ Jg* = 
0) is due to the global chiral symmetry(\[' e'''^"*!'). 

Different from classical mechanics, in the world of quantum 
mechanics, the chiral symmetry is broken In order to find 
what is happening, we need to analyze the vacuum structure: 
\vac\ Nl, Nr) = \vac; Nl) ® \vac; Nr), where 



To obtain the above results, we have regularized the diver- 
gent eigenvalues on the vacuum by ^-function regularization, 
for example, the gauge charge is regularized as follows: 



lim I > aj,a„ 
\nez 



+ 



nez 



(4.16) 



where A is an arbitrary constant with dimension of length 
which is necessary to make Ae^ dimensionless. This regular- 
ization respects gauge invariance because the energy of each 
level is a gauge invariant quantity. 

It can be shown that the gauge charge (Q) must vanish in 
order that the state is a physi cal st ate. We now have Nl = 
Nr. From the above equation (4.14), it can be seen that if A''^ 



and Nr are conserved, then, by varying A (the magnetic field), 
the chiral charge also changes. Therefore it is not a conserved 
quantity. We thus see that the vacuum is responsible for non- 
conservation of chira lity even though the dynamics is chirally 
invariant. From eq.(4.9) we see that the chiral current jf is 
equivalent to the electric current in the tubule axis direction. 



5 = f Js ix)dx = 



J^(x)dx = — J^. 



(4.17) 



Hence, in order to observe the anomaly, we should observe 
the electrical current in the torus. 

Due to the existence of the two Fermi points, the total 
current in the torus is given by the sum of two currents {Qs}^.^ 



and (Q s);.^ . We should care for the sign(±) 
charge (4.14). We define 



{Q5)k, ^ Nl+Nr-2N -1 + -N, 
{Q5)k,=NL+NR-2N-l-^N, 



in the chiral 

(4.18) 
(4.19) 



and treat them separately. 

It is clear from the above equations that there are two ori- 
gins of the usual current fiow in the torus. One is the Nl+Nr 
term which can be induced in thermal bath or by a sudden 
change of the magnetic field. On the other hand, magnetic 
field can change the quantum vacuum structure and lead to 
the anomaly. In order to avoid the unexpected changes of 
Nl{~ Nr), the magnetic field must be changed adiabatically 
at low temperature(< ^^). However, in the adiabatic pro- 
cess, when the strength of the magnetic field reaches the spe- 
cific points, then Nl{= Nr) also have to change. For simplic- 
ity, we set A'' = 1 and focus on the (Qs)^^. When increasing 
N starting from the point A = 0, Nl = 0, the energy is going 



\vac; Nl) = TT ai\0}, \vac; Nr) = TT bi\0}. (4.12) ^p as in eq(yj). At A = 3, the spectrum meets an another 



We define \vac; NL){\vac; Nr}) such that the levels with en- 
ergy lower than e]Vj^(— ejvj,-i) are filled and the others are 
empty. On this vacuum the charge expectation values and 
the energy become 



line of spectrum Nl = 1, see Fig.^. Therefore the circular 
current for fci in the ring 



J\ki) = 2{N - Nl) + ^ 



(4.20) 



follows the line shown in Fig.| 



(Q) =Nl- Nr, 

(Qs) = Al + A-fl - 2A- - 1 ± |a-, 



(Hf) 



27Th f {Qf + iQ^f 
L \ 4 



12 



(4.13) 
(4.14) 

(4.15) 
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FIG. 4. Adiabatic change of the fermionic energy of fci 
point. The energy value is labeled in the unit of 2i£: 




FIG. 5. Induced current in the twisted zigzag torus; 

The same analysis can be applied to the case of J^{k2) ~ 
— {Qb)f^^, the adiabatic change of energy and induced current 
are shown in Fig.n and Fig.[7l 



Hf 




FIG. 6. Adiabatic change of the fermionic energy of k2 
point. The energy value is labeled in the unit of 2i£ 




FIG. 7. Induced current in the twisted zigzag torus: 

Total current on the torus is given by a sum of two current, 
J = j\ki) + J^(k2). (4.21) 

Magnetic field dependence of the total current is shown in 
Fig^. Dotted line in Fig.^ indicates the current on the 
untwisted(Ai' — 0) torus. At each Fermi point, there are two 
spin degrees of freedom. Therefore the actual current is twice 
the J. This current for untwisted torus shows the same mag- 
netic field dependence to the persistent current in ref. [ p^ . 




FIG. 8. Magnetic field dependence of the total current on 
the torus: J = J^(fci) + J^(fc2) 



V. VACUUM STRUCTURE OF A CARBON 
NANOTORUS 

In this section we consider the vacuum structure of the total 
Hamiltonian: 



H 



He 



H^Hf + He, 



8tt 



\x - y\ 



dxdy. 



(5.1) 



In the previous section we have solved the kinetic part of the 
total Hamiltonian. We clarified its vacuum state and obtained 
the regularized eigenvalues of the physical quantities. We saw 
that the vacuum has the chiral charge so that it read the 
chiral anomaly. What we are interested in at this point is 
whether the previous results change or not by the inclusion of 
the Coulomb interaction in our analysis jl^ . Furthermore we 
hope to make clear the effects of the 'external' charge density 
on the chiral anomaly. Their understandings are the first 
step for studying more physically interesting situations such 
as impurity effects and junctions of a CNT and a metal or a 
superconductor. 

The Coulomb interaction consists of a product of the charge 
density. Therefore it is very convenient to rewrite the kinetic 
term Hf using the current operators. For this purpose, it is 
useful to introduce the left and right currents as follows: 



f{x) = Jl{x) + Jr{x), 



(5.2) 



N 



where Jl{x) = 'I'^(x)'I'l(x) and Jr{x) = *|j(a;)*H(x). We 
expand these currents by the Fourier modes, 
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Mx) = 



3l 



Jr{x) 



(5.3) 



where the Fourier components are the bosonic operators, 



(5.4) 



which satisfy the following commutation relations on the 
fermion Fock space, 



r .n / •m\'\ 
\jL,(J' ^ 



(5.5) 
(5.6) 



It is well-known that the following Hamiltonian has the same 
matrix element as the original fermion Hamiltonian: 



Hp 



L '\\ 4 12 



(5.7) 



where Vf '^s the Fermi velocity (v_f = \I3\). This Hamiltonian 



has a new term which is not shown in oq.(4.15). This term 
hcis vanishing value on the previous vacuum state because 



jl\vac\ Nl) = 0, jR\vac\ Nr) = 0, 



(5.8) 



for positive n. 

The Coulomb interaction can be rewritten using the 
bosonic current operators. 



He = 



— T 

n>0 



^^(n)((j2)t+jS)(jE + (jS)0- (5.9) 



Here, we introduce the Fourier component of the Coulomb 
potential. 



V{n) =2n dz 
Jo 



cos(27rnz) 



(5.10) 



where R is the circumference of the tubule. 

Some comments are in order. When we write the Coulomb 



interaction in eq.(5.1), it has an ultraviolet divergence in the 



limit oi X ^ y. Therefore we need to introduce some cutoff 
length. It is appropriate that we set it the length of a diameter 
of a tubule because we make the approximation that mixing of 
the different momentums in the compactified direction cannot 
happen. This explain the term in the denominator. 

Besides, the form of a torus is not a line but a ring on a plane 
so that we should use the direct length between x and y in 
the Coulomb interaction. This is the origin of the sin (tyz) 
term in the denominator. 

We combine the kinetic term and the Coulomb term as 
follows: 



n>0 



where 



(5.11) 



(5.12) 



and 



rr ^uTiVF r / .n\t -n . -n / ■n\\ 1 

Hr, = — J — {[JlYJl +Jr[Jr)' ~nj 



+- 



^^((jL)UiS)(jZ + (jS)t). (5.13) 



Here we introduce the fine structure constant a{— -^^)- 

We diagonarize the Hamiltonian Hn{n ^ 0) using the Bo- 
goliubov transformation: 



ri \ _ (^osht„ sinht„ ^| ( jl ^g^^^ 



(jS)V V™'^*" cosh t„ J y{jR 



where 



sinh 2tn 



1 e^T/(n) 



if 2nhvF eV(n) 
cosh 2t„ = — { + —Y- 



Er, = — r-^Wl + V(:n). 

L y TT Vf 



After some calculations we derive 



2-KhvF 



(5.15) 
(5.16) 

(5.17) 
n. (5.18) 



The energy E„ differs from 2ze£ilei The difference is due to 
the Coulomb interaction, see Figj^. 



L/R = 100 



En 



20 30 40 50 



FIG. 9. Energy spectrums are corrected by the Coulomb 
interaction. The energy val ue is labeled in the unit of '^^^^p . 
Here we set ^ = 10^ in eq.( |5.1C| ). 



The generators of the Bogoliubov transformation are given 



by 



Un = exp 



V2i 



-{{Jl) [Jr) - JlJr] 



Hence, we obtain the vacuum state as follows: 



\n>0 / 



\vac;NL,NR), 



(5.19) 



(5.20) 



where Nl = Nr. The previous vacuum state changes into 
new vacuum by the Coulomb interaction. We should estimate 
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the gauge charge and chiral charge on this new vacuum. Be- 
cause the operators(?7„) commute with the charge operators 
Q and Q5, the eigenvalues of these charges do not change. 



Q\vac) — 0, 



(5.21) 
(5.22) 



Hence, the chiral anomaly which we have considered in the 
previous section still exists when we include the Coulomb in- 
teraction into the analysis ^A^ . 



VI. EFFECTS OF EXTERNAL CHARGES 

We have so far considered the Coulomb interaction between 
the internal charges. However, more important problem may 
be the following: when we put external charges on a nan- 
otorus, how the vacuum structure change ? A fixed external 
charge may corresponds to an electrical contact or an impu- 
rity on a torus. In this section we put external charges on a 
torus and study the vacuum structure. Especially we analyze 
an effect of external charges on the chiral charge, the potential 
behavior between a pair of external charges and the 'charge 
screening' effect. The charge screening means that the inter- 
nal charge density is induced by the external charges so that 
the external charges are screened. 

We set two unit external charges on a torus. One has a unit 
charge e placed at xo and the other has an opposite charge 
— e at yo. 



Jex{x) = 5{x - Xo) - 5{x - yo) 



where 



Je 



(6.1) 



(6.2) 



We consider the Coulomb interaction with the following 
charge density which is a sum of the internal charges and 
the external charges: 



J°ix) + jUx) = 



(6.3) 



First we can easily show that the chiral charge on this new vac- 
uum does not change even in the presence of external charges. 



/ ■ tO I t5 / \\ tO \ (Os) 

{ya,c;Jex\J {x)\vac;J^^) = ^—. 



(6.7) 



So, the current by the chiral anomaly is not affected by the 
charged impurities. Second, we estimate the energy change 
due to the presence of the external charges: 

E{xo — yo) = {vac; j'ix\H{Jex)\vac\ J°^) — {vac\H\vac) 

E /3n / 2-KhvF \ ^ , .„ .„ 

n>0 



2-KhvF sr-^ 



a c 

TV Vp 



V{n) 



^ 1 + f^J^Vin) 



^^^2.n{xo-yo) ^^^^ 



where H{Jsx) is the Hamiltonian with the external 
charges( Jea;). Fig.^ shows the potential energy (eq.(3.5)) as 
a function of the distance between the two external charges. 



E(x) . 



L/R = 100 



FIG. 10. Distance dependence of the energy caused by two 
fixed external charges. The energy value is labeled in the unit 
of ^'^^"^ . The distance between two external charges is scaled 



by L. Here we take = 10 in eq.(p.l 

We see that the effect of the 'charge screening' on the ex- 
ternal charge can not be ignored in a metallic nanotorus. Be- 
cause the potential energy is now shown to be short-ranged. 
This means that some internal charges are influenced by the 
external charges, and external charges are screened. To con- 
firm this, we also compute the induced internal charges dis- 
tribution, 



{vac; J°^\f\x)\vac; J°^) = f(x; 



Xo 



fix;yo) 



(6.9) 



After some calculations we get 



2'Khvp ^ f3n /2nhvF 



L 



{jex)*jex, (6.4) 



where 7„ — sinh2t„(coshf„ — sinht„) and /3„ — ■£;^V{'n)- It 
is easy to find conditions of the vacuum in the presence of the 
external charges, 

(ii + Jnjex) Ivac; jI^) = 0, (6.5) 
{fR + 1riUex)*)\vac;J°^)=0, n > 0. (6.6) 



where 



f{x;xo) = -y 



a c_ 

TV Vp 



V{n) 



2Tin(x — Xo) 



n>0 



:V{n) 



(6.10) 



This function is displayed in Fig.pj]. 

We should stress here that t he ab ove analysis on the charg- 
ing energy ( ^.^ ) and screening ( |6.10| ) are not complete because 
the charge density {j'^(x)) in the Coulomb interaction consist 
of one massless fermion in our analysis. A CNT have four 
independent fermions due to the two Fermi points( 2.18 ) and 
spin degrees of freedom. However, straight forward extension 
of our anajvsis shows that the conclusions in this paper hardly 
changes 
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L/(x;0) 







L/R = 100 



[x/L] 

FIG. 11. The internal charge density modulation around a 
unit charge e placed at origin: /(x;0). The charge density 
value and the distance from the charge are labeled in t he un it 
of L and j- respectively. Here we take — 10^ in eq.(5.1C). 



VII. CONCLUSION AND DISCUSSION 

In this paper we analyzed the carbon nanotorus and dis- 
cussed the quantum mechanical vacuum structure of a metal- 
lic torus. The significance of our present work may be put as 
follows. 

We used the quantum field theory to analyze the vacuum 
structure of a metallic nanotorus. We pointed out that the 
chiral anomaly in 1-1-1 dimensions should be observed in the 
form of specific magnetic field dependence to the current. 
This current is the same as the persistent current. It is cer- 
tain that the persistent current can be understood in the light 
of the chiral anomaly in 1-1-1 dimensions. We also clarified 
the vacuum state including the Coulomb interaction and dis- 
cussed the effect of 'charge screening' on the external charges. 
It was found that the chiral anomaly is not affected by the 
charged impurities and the charge screening effect actually 
occurs. 

Although we have analyzed a nanotorus in the present pa- 
per, we can regard the torus as a tube locally provided the 
circumferential length of the torus is large. Hence, we may 
apply the results obtained in the analysis of the torus to the 
study of a carbon nanotube. To take an example, it is rea- 
sonable to suppose that the screening effect is also present in 
a tube like in a nanotorus. 

What we do not consider in this paper are effects of finite 
temperature( > '^'^^^f ) on chiral anomaly and screening of 
external charges. In finite temperature, in addition to the 
electron field, phonon field would come into play |^^. This 
phonon field may be vital to understand the electrical behav- 
ior of the carbon nanotube at finite temperature. We would 
like to make a quantitative analysis of finite temperature ef- 
fects in a future report. 
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